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Seesaw Scale in the Minimal Renormalizable SO (10) Grand Unification 
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Simple 50(10) Higgs models with the adjoint representation triggering the grand-unified symme- 
try breaking, discarded a long ago due to inherent tree-level tachyonic instabilities in the physically 
interesting scenarios, have been recently brought back to life by quantum effects. In this work we 
focus on the variant with 45 h © 126 h in the Higgs sector and show that there are several regions in 
the parameter space of this model that can support stable unifying configurations with the B — L 
breaking scale as high as 10 14 GeV, well above the previous generic estimates based on the minimal 
survival hypothesis. This admits for a renormalizable implementation of the canonical seesaw and 
makes the simplest potentially realistic scenario of this kind a good candidate for a minimal SO (10) 
grand unification. Last, but not least, this setting is likely to be extensively testable at future 
large-volume facilities such as Hyper-Kamiokande. 



PACS numbers: 12.10.-g, 12.60.Jv, 12.15.Ff 



I. INTRODUCTION 

For the last about thirty years, the simplest non- 
supersymmetric 50(10) gauge models with 45// © 16// 
or 45// © 126// in the Higgs sector have been widely 
considered uninteresting for any realistic unified model 
building. This was namely due to the tachyonic insta- 
bilities in their tree-level spectra popping up in all set- 
tings compatible with the basic gauge unification con- 
straints [l|-|5[ which, in non-SUSY settings, generically 
favour intermediate-energy thresholds. However, as it 
was shown recently in |,|, such instabilities are just 
artefacts of the tree-level approximation. Hence, tech- 
nically, quantum effects bring this class of models back 
from oblivion. 

On the other hand, dedicated renormalization group 
studies such as [8l-4ll1] reveal that a successful unifica- 
tion in this class of models typically requires the B — L 
breaking scale below 10 12 GeV for the 45//© 16// variant 
and below 10 10 GeV in the 45// © 126// case. Such val- 
ues, however, are disfavoured by the neutrino oscillation 
and cosmology data: i) In the former case, (16//) breaks 
the B — L symmetry by one unit and, thus, the seesaw 
requires a pair of (16//) insertions. This can be mini- 
mally implemented at the renormalizable level by e.g. a 
variant of the Witten's radiative mechanism (12h14| or, 
giving up renormalizability, by a d = 5 operator. In ei- 
ther case the "effective" A(B — L) = 2 seesaw scale is 
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further suppressed with respect to the B — L breaking 
scale and the light neutrino masses are typically over- 
shoot by many orders of magnitude. Moreover, the non- 
rcnormalizable nature of the seesaw in the d = 5 case 
hinders the general predictivity of this model, ii) With 
126// at play, the B — L symmetry is broken by two units 
so the right-handed neutrinos receive their masses at the 
tree level via the renormalizable 16/?16f126^ Yukawa 
interaction [HI, [lf3]- The upper limit on (126//) quoted 
above then again pushes the absolute scale of the light 
neutrino masses much above the current limits. 

Though unpleasant, this, however, does not constitute 
a fundamental blow to the minimal non-SUSY 5*0(10) 
as an extensive multi-parameter fine-tuning in the see- 
saw formula can still bring the light neutrino masses 
down to the desired sub-eV domain. In this respect, 
the situation is very different from that of the mini- 
mal supersymmetric (SUSY) 5*0(10) grand unified the- 
ory (GUT) p7H25| where the neutrino masses are typ- 
ically undershot; indeed, the rigidity of the Higgs po- 
tential in minimal SUSY Higgs models enforces a pop- 
ulation pseudo-Goldstonc bosons well below the GUT 
scale (M G ) whenever the 50(10) SM breaking 
is not essentially one-step [2(1 H3] , hence disturbing the 
nearly ideal unification within the Minimal Supersym- 
metric Standard Model (MSSM); for further information 
see, e.g., [28| and references therein. 

In the same spirit, one should keep in mind that the key 
upper bounds on the B — L scale identified in [8MTT1| are 
derived under the strong assumption of the minimal sur- 
vival hypothesis [29j, i.e., that a minimal set of needed 
intermediate thresholds cluster exactly at the relevant 
symmetry breaking scale. This, of course, does not need 
to be the case in general and as little as a single unex- 
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pected multiplet in the bulk can open a room for B — L 
scales much above the naive expectation, thus rendering 
the gauge coupling unification compatible with the neu- 
trino data for a reasonable price. In this respect, the 
non-SUSY models with higher-dimensional Higgs repre- 
sentations (such as 45 h © 126//) featuring a number of 
free parameters in the Higgs potential provide a lot of 
room for such a serendipity. Moreover, given the rcnor- 
malizable nature of the seesaw in the 45// © 126// case, 
the Yukawa sector of this kind of models is strongly con- 
strained, which further opens the door for their near fu- 
ture testability. 

In this study we focus on the possible role of acci- 
dental thresholds in the desert of the minimal SO(10) 
GUTs based on the 45// © 126// Higgs sector. In par- 
ticular, we calculate the tree-level spectrum of the mini- 
mal Higgs model and the leading universal radiative cor- 
rection to the relevant Higgs masses and ask ourselves 
i) whether states with accidentally small masses can pop 
up in some regions of the parametric space without desta- 
bilising the scalar potential and it) whether the corre- 
sponding threshold effects can lift the seesaw scale to the 
desired ballpark of 10 13 - 14 GcV. 

The work is organized as follows: In Section [TT] we de- 
fine the 45// © 126 ff 50(10) Higgs model of interest and 
calculate its tree-level spectrurrQ which reveals the ex- 
pected tachyonic instabilities except for the phenomeno- 
logically questionable SU (5)-like descents. In analogy to 
the canonical example elaborated on in [(| we argue that 
radiative corrections alleviate the issue and that stable 
and potentially realistic SM vacua are accessible. To ex- 
emplify that, we calculate the leading 5O(10)-invariant 
radiative correction as a minimal scalar-spectrum regula- 
tor. In Scction lTHl we study the possible effects of various 
multiplets - if they happen to live in the "GUT desert" 
- on the actual location of the B — L scale. We iden- 
tify two specific simple and consistent settings in which 
all current phenomcnological constraints from the proton 
decay searches and big-bang nucleosynthesis are compat- 
ible with the latest limits on the absolute neutrino mass 
scale. A simple numerical scan over the parametric space 
reveals extended domains supporting these solutions. Re- 
markably enough, in both cases the extra threshold is 
pinned to a relatively narrow mass window which, in 
turn, yields a rather specific prediction for the position of 
the GUT scale and, hence, the d = 6 proton decay rate, 
well within the reach of the future large volume facilities 
such as Hyper-Kamiokande (Hyper-K) (32[. 



With all this at hand, in Section IIVI we make a case 
for a new potentially realistic minimal renormalizable 
50(10) GUT based on the 45// © 126// © 10// Higgs sec- 
tor. We comment in brief on the prospects and strategies 
of a future more detailed scrutiny of the scheme, paying 
particular attention to the Yukawa sector fits and the 
ultimate calculation of the proton decay branching ra- 
tios in the fully consistent settings. Then we conclude. 
Technical aspects of the Higgs and gauge-boson spectrum 
calculation are deferred to a set of Appendices. 



II. THE 45-126 HIGGS MODEL 
A. The tree-level scalar potential 

The most general renormalizable scalar potential that 
can be written with 45// and 126// at hand reads 

V = U45 + U 126 + Knix , (1) 

where 

V i5 = ~y (#)o + y(#)o(#)o + y(#M#) 2 , (2) 

Ul 26 =-|W)o (3) 

+ ^.(££*) (££*)o + ^(EE*) 2 (EE*) 2 

+ (3!)W (SS * )4(SS * )4 + W (SS * )4 ' (EE * )4 ' 
+ (^(SS) 2 (ES) 2 + ^(S*S*) 2 (S*£*) 2 , 

U mix = J(0) 2 (SS*) 2 + _£L(^) (EE*) (4) 

+ ^|j(#)4(SS*)4 + §(#) 4 <(££> 

+ ^(#) 2 (££) 2 + |[(#) 2 (£*£*) 2 . 

Here we have used the symbols <p and £ for the compo- 
nents of 45// and 126//, respectively. The detailed break- 
down of all the contractions (with the subscripts denoting 
the number of open indices in the relevant brackets) is 
given in Appendix [A] Notice that all the couplings are 
real but r\i and 72. 



B. The symmetry breaking patterns 



1 Here the non-SUSY nature of the model is central - the SM- 
vacuum manifold of the minimal SUSY GUT, as complicated 
as it naively looks, is in reality very simple; indeed, it is 
parametrized by a single complex parameter [20j . 

2 Though there exist detailed studies of the vacuum of the 54// © 
126// Higgs model, cf. |30ll3l|. a similar analysis for the setting 
with 45// instead of 54//, to our best knowledge, has never been 
done. 



There are in general three SM singlets in the re- 
ducible 45// © 126// representation of 5O(10). Using 
BL = (B — L)/2 and labelling the field components 
with respect to the 3 c 2 l 2 r 1bl (i.e., 5£/(3) c ©5t/(2) L © 
SU(2)n<g) U(1)bl) algebra, the SM singlets reside in the 
(1, 1, 1, 0) and (1, 1, 3, 0) sub-multiplcts of 45// and in the 
(1, 1, 3, +1) component of 126//. In what follows we shall 
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denote 

((1, 1, 1, 0)) = uj BL , ((1, 1, 3, 0)) = u R , ((1, 1, 3, +1)) = a, 

(5) 

where ujbl.r are real and a can be made real by a 
phase redefinition of the 126//. Different VEV configura- 
tions trigger the spontaneous breakdown of the 5*0(10) 
symmetry into several qualitatively distinct subgroups. 
Namely, for a = one finds (in an obvious notation) 



of pseudo-Goldstone bosons (cf. also comments in Ap- 
pendix |B4]) entertaining very simple mass formulae: 



LU R = 0, UJ B L ^ 

UJr ^ 0, UJ B L = 

UJR ^ 0, UJRL ^ 

UR = -UBL 7^ 

ur = ujbl ^ : 



3 C 2l 2r Ibl ) 

3 C 2l 1r Ibl , 
flipped 5' lz' , 
standard 5 lz , 



(6) 



with 5 lz and 5' lz' standing for the two inequivalent 
embeddings of the SM hypercharge operator Y into 
SU(5) (g> U(l) C 50(10) usually called the "standard" 
and the "flipped" SU(5) scenarios [H,[33], respectively. 

(3) 

In the standard case, Y = T R +T BL belongs to the 
SU(b) algebra and the orthogonal Cartan generator Z 

is given by Z = -4T { R 3) + 6T BL - In the flipped (5'1 Z ,) 
case, the right-handed isospin assignment of quarks and 
leptons is turned over so that the flipped hypercharge 
generator reads Y' = —To 3 ' +Tbl- Accordingly, the ad- 
ditional U(l) z > generator reads Z' = AT R 3) + 6T BL (for 
further details see ,e.g., Ref. @). 

For (7^0 all the intermediate gauge symmetries flBJ) 
are spontaneously broken down to the SM group, with 
the exception of the last case which maintains the SU (5) 
subgroup unbroken and, hence, will not be considered 
here. The decomposition of the 45// and 126// repre- 
sentations with respect to the all relevant intermediate 
symmetries ^ is detailed in Tables HI and IPTl 



C. The tree-level scalar spectrum 

Adopting the convention in which the mass term in the 
Lagrangian is written as \i\j t M 2 ^, where ip = (</>, £*, E) 
is a 297-dimensional vector, the scalar spectrum is ob- 
tained readily by evaluating the relevant functional scalar 
mass matrix of the schematic form 



M 2 (0,E*,E) 




(7) 



on the SM vacuum. The subscripts here denote the 
derivatives of the scalar potential with respect to a spe- 
cific set of fields. Subsequently, this matrix is brought to 
a block-diagonal form by a subsequent unitary transfor- 
mation into the SM basis. 

The complete tree-level spectrum is given in Ap- 
pendix [Bj There are several features that can be seen 
readily: i) as anticipated in Q there is again a pair 



M 2 (l,3,0) 
M 2 (8,l,0) 



-2a 2 (uj B L - ur)(uj B l + 2ujr) , (8) 
-2a2(ujR~ ujbl)(ujr + 2uj B l) ■ (9) 



These multiplets develop tachyonic masses whenever 
ubl/wr is outside the [— 2, — |] interval. Hence, as such, 
the tree-level Higgs spectrum is clearly unable to support 
the physically interesting breaking patterns with cither 
ubl "C ujr or a;^ <C ujbl, thus avoiding the intermedi- 
ate flipped SU(5)' ® U(l)z' stage. SU(5) intermediate 
stages, ujbl,ur <C <t, are disfavoured as welfl. 

ii) In this respect, it is worth looking at formulae (|8|) 
and (EH) in more detail. For instance, as in the 45jj © 16# 
case Q there are no contributions there from the B — L 
-breaking VEV a although the number of available con- 
tractions of the type (0 2 )(SI]'*^) is larger here. This can 
be understood as follows: regardless of how the indices 
of the EE 1 -*) bilinears are contracted, the resulting ten- 
sor never breaks the SU(5) symmetry. Since both E's 
couple to both 0's in the same manner one can always 
view the contraction with the pair of adjoints (0's) as 
a quadratic covariant-derivative-like term for the fields 
with the SM gluon and SU (2) ^-gauge quantum numbers. 
These fields, however, remain massless at the SU (5) level. 

However, as shown in Q, this is no longer the case 
at the quantum level where all the tree-level forbidden 
couplings do indeed enter the relevant mass formulae and 
thus open the room for the physically interesting settings 
with ujrl very different from ujr. 



D. Leading one-loop corrections 

Unlike in the 45 h © 16// case the full-fledged effective 
potential (EP) calculation of the one-loop scalar spec- 
trum in the 45# © 126 h model is very difficult due to the 
enormous complexity of the contractions involving 126# 
so we shall not attempt it here. However, the radia- 
tive corrections are really important only for the pseudo- 
Goldstone bosons associated to accidental global symme- 
tries, cf. 0] and Appendix IB 4[ thus, one can get a good 
grip on the one-loop spectrum even without the full EP 
analysis. Moreover, some of the results obtained for the 
45/fffil6# setting in @ can be readily adopted to the cur- 
rent case; in particular, the one-loop gauge-induced cor- 
rections to the masses of the scalars residing solely in 45// 



3 Recently, there were several attempts to reconcile the simplest 
non-SUSY SU (5) scenarios with the gauge unification by means 



of intermediate-scale thresholds, see, e.g., I35H41H and references 
therein. Though possible in principle, we shall not consider the 
SU(5) option here because there is virtually no room for such 
an intermediate stage below the 5O(10)-breaking scale in the 
model of our interest. Moreover, these settings generically rely 
on several multiplets pushed into the desert. 
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TABLE I. Decomposition of the adjoint representation 45 with respect to the various SO(10) subgroups. The definitions and normalization 
of the abelian charges are given in the text. 
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(6,1,1) 


(6,1,0) 
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(3,1,0,+i) 


(3,1,+!) 


(5, +2) 


(5, +2) 
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(3,1,-!) 


(45,-2) 
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(TO, 1,-1) 
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(1,1,-1,+1) 


(1,1,0) 


(1,+10) 


(50, +2) 




(10,1,0) 




(1,1,0,+1) 


(1,1, +1) 


(10, +6) 


(10, +6) 




(T0,l,+1) 




(1,1,+1,+1) 


(1,1, +2) 


(50, +2) 


(1,+10) 






(3,1,3,+|) 


(3,1,-1,+D 


(3,1,-1) 


(10, +6) 


(45, -2) 
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(3,1,+!) 
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(50, +2) 








(3,1,+1,+|) 


(3,1,+i) 
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(10, +6) 






(6,1,3,-i) 


(6,1,-1,-!) 


(6,1,-!) 


(50, +2) 


(15, -6) 








(6,1,0,-i) 


(6,1,-!) 


(45,-2) 


(45,-2) 








(6,1, +1,-|) 


(6,1,+!) 


(15, -6) 


(50, +2) 


(15,2,2) 


(15, 2, -i) 


(1,2,2,0) 


(1,2,-1,0) 


(1,2,-!) 


(5, +2) 


(45, -2) 




(15,2,+i) 




(l,2,+i,0) 


(1,2,+!) 


(45,-2) 


(5, +2) 






(3,2,2,-f) 


(3, 2, — j , — 3) 


(3,2,-1) 


(45,-2) 


(15,-6) 








(3, 2, +1, — |) 


(3,2,-!) 


(15, -6) 


(45, -2) 






(3, 2, 2,+|) 


(3,2,+i+f) 


(3,2,+D 


(50, +2) 


(10, +6) 








(3, 2, — |, +|) 


(3,2,+!) 


(10, +6) 


(50, +2) 






(8,2,2,0) 


(8,2,-1,0) 


(8,2,-!) 


(50, +2) 


(45, -2) 








(8, 2, +1,0) 


(8,2,+!) 


(45,-2) 


(50, +2) 



TABLE II. Same as in Table U for the 126 representation. 



(such as (8,1,0) and (1,3,0)) are identical to those ob- 
tained in Q, cf. formulae (D1)-(D2) therein. This, how- 
ever, is not the case for the contribution of scalars which 
span over the components of 126 h- Similarly, the one- 
loop scalar-induced contributions to the tree level scalar 
masses should be calculated from scratch. Needless to 
say, this is a formidable task if it is to be performed in 



full generality. 

Thus, in what follows, we shall focus only on the most 
universal scalar one-loop correction, namely, the leading 
non- logarithmic 5'O(10)-invariant r 2 -proportional term 
which, as we argue, can be fully accounted for by a sim- 
ple diagrammatic calculation. Since it yields a positive 
correction to all the scalar masses, it should already be 
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enough to regularize the salient tachyonic instabilities of 
the tree-level scalar spectrum and, perhaps, open new 
regions in the parametric space where stable unifying 
configurations with phcnomcnologically favourable inter- 
mediate scales could be supported. Moreover, since also 
the other leading non- logarithmic corrections (i.e., those 
coming from the gauge and the remaining scalar loops) 
are typically positive, including just the SO( 10) invariant 
piece can be viewed as a minimalistic attempt to stabilize 
the tachyons. In view of this, a detailed calculation of all 
one-loop corrections to the scalar spectrum in this frame- 
work is not even necessary and will be left to a dedicated 
future study. 

Since the leading scalar-loop induced non-logarithmic 
corrections in the scalar sector come from tadpoles [42j . 
it is easy to see that the only source of a r 2 -proportional 
non-log term is associated to the renormalization of 
the stationarity conditions. Diagrammatically, it corre- 
sponds to a special cluster of one-loop graphs contribut- 
ing to the one-point function of 45# of the kind 




-o g 



;---o 



r ;---o+. 



126' to 



(10) 

Given the S'O(IO) structure of the relevant r-vertex, 
namely, 



IT 



U45-I26 3 ^^ij^klmni^klmnj 



(11) 



the universal mass shift due to this class of graphs reads^l 

35r 2 



AM? , 2 = 

1-loop-r* 327r 2 



logs, 



(12) 



where the symbol "logs" denotes all the logarithmic cor- 
rections that arc minimized at the GUT scale. 



III. UNIFICATION IN THE 45-126 MODEL 

With this information at hand, in this section we can 
finally address the question of our main interest, namely, 
whether accidentally light scalar multiplets in the SM 
desert could possibly open the door to a consistent gauge 
unification with a B — L scale well above the unpleasant 
upper limit of about 10 10 GeV obtained in [ll| under the 
assumption of minimal survival. 

Since the scalar masses are expressed as functions 
of the microscopic parameters entering the scalar po- 
tential ([I}, pushing a specific multiplet into the desert 



amounts to imposing an extra algebraic constraint on 
the parameter space of the model, i.e., it cuts out a re- 
gion close to the relevant zero-mass hyper-surface. The 
rest of the spectrum then must be evaluated around this 
hyper-surface which, however, brings in a high level of 
non-linearity. Thus, in what follows, we shall mainly 
stick to numerical methods to simulate the heavy scalar 
and vector-boson spectra in order to single out the re- 
gions of the parametric space that can support viable 
gauge unification patterns. 



A. Consistency 

Besides gauge unification, there are other basic aspects 
of an overall consistency of potentially realistic settings 
that will be of our concern here, namely, the stability of 
the physical vacuum (i.e., the absence of tachyons) and 
the position of the unification scal^3 which governs the 
d = 6 proton decay. Moreover, with potentially very light 
coloured states in the desert, d > 6 proton decay as well 
as possible BBN issues should be also considered. 



1. Vacuum stability 

As stated before, from now on we shall stick to 
the "minimally regularized" form of the scalar spec- 
trum, i.e., we shall use the tree- level formulae of Ap- 
pendix [22] augmented with the leading S'O(10)-invariant 
non- logarithmic one-loop correction (|12[) . 

For each physical point in the parametric space, the 
mass-squares of all propagating scalars should be posi- 
tive. This, as we shall sec, is indeed a very restrictive 
constraint which already disqualifies some of the poten- 
tially interesting multiplets, see Sect. IIIIB 11 It is per- 
haps worth mentioning that with one such a stable vac- 
uum at hand one can generate a continuum of other sta- 
ble vacua by rescaling all the dimensionful parameters 
entering the mass formulae by a common factor. This in- 
variancc will be later on used for a simple optimization of 
the one-loop unification patterns, cf. Sect. MI B 21 More- 
over, it is easy to understand that, as long as only the 
scalar mass-squares are concerned, further degeneracies 
in the parametric space of the model can be identified; 
among these perhaps the most prominent is the absence 
of the phase of 72 from the tree-level mass formulae and 
the irrelevance of the overall sign of the mass parameters 
at play (i.e., all that matters are just relative signs). 



4 Let us just mention that the same technique applied to the 
45j^ ffi 16h case yields the uniform mass shift of T 2 /4-7r 2 which 
is indeed consistent with the results of the effective potential 
calculation [||. 



5 The GUT scale is conventionally defined as the mass scale of 
the gauge bosons associated to the breakdown of SO(fO) to ei- 
ther Zc'Zl'ZrXbl ° r ^C^L^R intermediate symmetries, i.e., those 
transforming for instance as (3, 2, — ^) © (3, 2, +§) under the SM 
gauge group. 
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2. Proton lifetime limits 



3. BBN and the lifetime of light coloured BSM multiplets 



a. d — 6 proton decay: We shall impose the latest 
(2011) Super-Kamiokande (SK) limit on the proton life- 
time (for the e + 7r° channel) |43j : 

r{p -)■ e+7r°)sK,20ii > 8.2 x 10 33 years , (13) 

and, whenever appropriate, comment on the changes in 
the results for a couple of assumed future sensitivity lim- 
its, namely those quoted in [32| that Hyper-Kamiokande 
(HK) should reach by 2025 and 2040, respectively: 

r(p^ e + 7r°) H K,2025 > 9 x 10 34 years, 



r(p^ e + 7r°) H K,2040 > 2 x 10 30 years 



35 , 



(14) 
(15) 



These translate to the following (raw) formula for the 
compatibility regions in the Mg — ctQ 1 plane: 



45 



10 



2(n G -15) 



> 11.8, 39.0, 58.1, 



(16) 



where tig = log 10 (MG/GeV) and the three values on the 
right-hand-side correspond to the three lifetime limits in 
Eqs. (|13p - (fT5l) . respectively. In the relevant figures (cf. 
FIGs [TJS1 and FIGs , the regions of the parametric 
space where the three constraints (fTo) are fulfilled will 
be, consecutively, denoted by light-gray, dark-gray and a 
black color. 

One should also check that lowering a specific multiplet 
into the GUT desert does not bring any of the proton- 
dangerous coloured scalar triplets too much below some 
10 14 GeV; although the detailed structure of the scalar 
d = 6 proton decay amplitude is typically suppressed by 
small Yukawa couplings, this is not always the case and 
a coloured triplet well below this limit can be dangerous. 
Since we do not consider the details of the Yukawa sec- 
tor here, we shall adopt a conservative limit like the one 
quoted above. Remarkably enough, this constraint turns 
out to be rather weak and in a vast majority of the cases 
where (fT6|) are obeyed the scalar triplets are safe. 

b. d > 6 proton decay: Under the "big desert" hy- 
pothesis the d = 6 proton deca y ope rators conserve B — L 
up to Mw/Mg corrections [44ll45l pl. However this picture 
does not need to hold anymore if we consider new struc- 
tures at intermediate scales well below Mg and d > 6 pro- 
ton decaying operators (such as those conserving B + L 
at the d = 7 level, cf. [4(| |47j) should be inspected. 
A "canonical" example here is the situation when the 
(3,2,+i) scalar approaches the weak scale; the relevant 
B+L conserving proton decay amplitude can then easily 
clash with the experimental limits [H[ . 



6 In the SO(10) models these operators are usually induced by 
the scalar triplets transforming as (3, 1, — 3)® (3, 1, +5) and the 
(3,2,-|)®(5,2,+|)e(3,2,+|)e(3,2,-|) gauge bosons. 

7 In the current S'O(IO) model the relevant effective operator is 
traced back to the 126|j quartic coupling and the 16pl6p 126^ 
Yukawa interaction. 



Light colored thresholds can be also troublesome for 
the Big Bang Nucleosynthesis (BBN). This has to do with 
the requirement that any colored state other than the SM 
fields must decay with a lifetime shorter than about 1 sec- 
ond, in order to preserve the classical predictions of the 
light elements' abundances |43j . From this perspective, 
renormalizable Yukawa couplings of such light scalars to 
the SM matter fields are welcome as the relevant decay 
widths arc typically large enough to be safe. 



B. Running with extra thresholds in the desert 

For the sake of simplicity, in what follows we shall 
entirely stick to the case with a single extra SM sub- 
multiplet of 45# © 126j/ in the desert. This not only 
lowers the number of fine-tunings to the minimum, but 
also admits for a systematic classification of the possible 
threshold effects. 



Identifying the most suitable thresholds 



i) The stability requirements of Sect. IIII A II disfavour 
a light (3, 3, — |) multiplet as there are no suitable stable 
vacua supporting this configuration even if the leading 
universal one-loop correction (fTTj) is taken into accounfl 

ii) There is a good reason to disfavour all multiplets 
whose effect on the hypercharge coupling evolution is 
much larger than the effect on the SU{2)l coupling: Re- 
call that the upper limit on B — L emerges from the 
need to delay the "premature" SM unification of the 
U{1)y and SU(2)l couplings by lowering enough the 
B — L scale. An extra state in the desert which would 
act against this rule of thumb would further strengthen 
the demands imposed on the B — L scale, thus further 
lowering the relevant upper bound. On the other hand, 
such states are almost never brought down alone as the 
relevant fine-tuning lowers also the states occupying the 
same larger-symmetry multiplets to the respective sym- 
metry breaking scale; however, such intermediate scales 
in the settings of our interest should not be far from Mg 
so the effects of these extra components are typically 
sub-leading. Hence, multiplets like (1,1, +1), (1,1, +2), 
(3, 1, -f ), (3, 1, +|), (6, 1, -|), (6, 1, - i), (6, 1, +|) and 
(3, 2, +|) are not fit for our purposes. On the same foot- 
ing, the individual effect of an additional (1,2,+^) is too 
weak to make much difference even if it is pushed down 
to the electroweak scale. 



Strictly speaking, this is not entirely decisive as the other one- 
loop corrections we are not considering here may open more room 
for such a setting. 
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Hi) We discard also the (1,3,-1) component of 126# 
because it is the type-II seesaw triplet - indeed, a very 
light triplet would require an extra fine-tuning of the 
effective S'C/(2)i-triplct-doublet-doublet coupling other- 
wise the absolute neutrino mass scale would be overshot 
by many orders of magnitude. 

Thus, from now on we shall entirely focus on the possi- 
ble threshold effects due to the remaining SM multiplets 
pushed into the GUT desert, namely, the (1,3,0) sub- 
multiplet of 45#, a pair of (3,2,+g) mixed multiplets, 
the (6, 3, +§) of 126 H and the pair of (8, 2, +5) in 126 H ■ 



2. Technical details of the RGE analysis 

On the technical side, we shall always work in the effec- 
tive SM picture where all the beyond-SM scalar and vec- 
tor bosons are classified by their 3 c 2l1y quantum num- 
bers; hence, conveniently, we will be always using the 
three SM "effective" couplings irrespective of the actual 
number of simple gauge factors that can be identified at 
any given energy scale. Needless to say, this is a mere 
convention provided that the matching to the full theory 
(especially at higher orders) is performed consistently. 

Given the tachyonic nature of the tree-level spectrum 
in the settings of our interest, a pure one-loop RGE anal- 
ysis is meaningless; in principle, the simplest fully con- 
sistent approach would be, of course, a two-loop run- 
ning based on a complete one-loop information about the 
scalar and gauge spectra. 

This, however, is extremely demanding in full gener- 
ality because even the very analytic minimization of the 
relevant one-loop effective potential in the 45# © 126h 
case is virtually intractable (note that, in this respect, 
the qualitative difference between the 45# © 16h and 
45# © 126r cases is paramount). 

Thus, we shall rather perform a qualitative one-loop 
RGE analysis based on the "minimally regularised" 
scalar spectrum, see Sections llI Dl and lHI A ll which. how- 
ever, should^] account for all the salient features of the 
fully consistent picture. In other words, we work in the 
approximation in which the full one-loop approach to the 
gauge coupling evolution is refined by the key two-loop 
effects. 

Technically, the calculations are performed in three 
stages along the following lines: First, we randomly scan 
over the parametric space of the model assuming the de- 
sired multiplet to be close to the electroweak scale and 
calculate the scalar spectrum for each such a point. For 



9 This expectation is based on the simple fact (see for instance [f| ) 
that the tachyons, which can be identified with pscudo-Goldstone 
modes of accidental global symmetries of the scalar potential re- 
stored in certain corners of the parametric space, are the only 
states whose mass-squares are really prone to radiative correc- 
tions. 



those points for which the vacuum turns out to be sta- 
ble, we adjust the overall scale of the dimensionful para- 
meters ujr, lobl, cr and r and the position of the light 
threshold in such a way that a consistent unification is 
obtained. Note that, in many cases, this can be done 
even analytically - at one loop, both such changes in- 
flict essentially linear shifts in the values of the three 
gauge couplings evaluated above the highest threshold 1 "*! 
so the optimization of the unification pattern amounts 
to a solution of a linear system. Finally, we check the 
full consistency of the resulting pattern with the proton 
decay and BBN limits specified in Section IIII Al and see 
whether the threshold effects can lift the B — L scale into 
the seesaw-favoured domain of lO 13 "^ 14 GeV. 

Given that, one can identify the following main sources 
of uncertainties plaguing the precision of the derived 
electroweak-scale values of the gauge couplings: i) Stick- 
ing to the one-loop beta-functions we commit an error 
of the size of a typical two-loop effectM. Assuming the 
usual size of such an uncertainty as observed, e.g., in [ll[ 
one can expect a reduction of the tree-level prediction of 
Mq by roughly a factor of two. ii) We do not re-input 
the derived values of the gauge couplings back into the 
gauge-boson mass-formulae and reiterate the code; for 
the sake of simplicity, we rather use a "typical" gauge 
coupling corresponding to a^} around 40 for all heavy 
gauge boson masses. The error due to this is of the or- 
der of loglsLc/sLumcd]/ 1671 " 2 which, however, is entirely 
negligible, in) We treat the fine-tuning in the (1, 2, +i) 
sector in a simplified manner: since one eigenstate of 
the (1,2, + 1) system is implicitly assumed to be fine- 
tuned to the electroweak scale to act as the usual SM 
Higgs boson only the heavier eigenstate of the doublet 
mass matrix should be included in the heavy-spectrum 
analysis. On the other hand, it does not make sense to 
perform such a fine-tuning with just 45# © 126h at play 
as it would, artificially, bring in an extra constraint on 
the parametric space which would be, however, absent 
in any realistic model including e.g. an extra 10h in the 
Higgs sector. For the sake of this qualitative analysis, we 
decided to resolve this dichotomy by mimicking the effect 
of the (unidentified) heavy doublet by averaging] over 
the effects of the two massive eigenstates of the doublet 



Concerning the latter, i.e., the shift in the mass of the acciden- 
tally light multiplet, as long as it is well below the next-to-lightcst 
and other heavier thresholds it can slide essentially freely with- 
out affecting the heavy part of the spectrum at all because all 
such configurations fall into a very small patch of the parametric 
space. This, however, does not need to be always the case and 
one should be more careful here, cf. FIG. [4] 

11 More precisely, here we refer to the typical scale of the two-loop 
corrections in "regular" settings, i.e., those in which one does not 
encounter tachyonic instabilities and a one-loop analysis would 
be perfectly self-consistent. 

12 Note that this way one implicitly retains the desired slope (6; = 
—37/3) for all three effective SM couplings above the scale of the 
heaviest component of the scalar spectrum. 
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mass matrix (|B9[) . Note also that the effect of a possible 
extra 10# in the full-fledged models (like those discussed 
later in Sect. II Vj) is expected to be small because the ex- 
tra degrees of freedom would typically cluster around the 
GUT scale and, amounting to a full irreducible 50(10) 
representation, they would affect the GUT-scale position 
only marginally. Last, but not least, an extra 10 h at play 
does not contain a new candidate for a suitable low-scale 
threshold so, in this respect, the classification given in 
Sect. IIIIBTl is not affected. 
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C. Results 

Let us begin with a short comment on the first two 
options identified in Sect. IIII B 11 namely, a light (1, 3, 0) 
and/or a light component of the (3, 2, +g) scalar pair. Al- 
though in both cases one can find regions of the paramet- 
ric space supporting such light multiplets in the desert, 
the predicted position of the GUT scale is always at least 
an order of magnitude below the current Super-K limit 
(cf. Sect. IIII A 2[) and, hence, the d — 6 proton decay con- 
straints are always badly violated. Thus, none of these 
two cases turns out to be interesting. 

However, as we argue below, in the remaining cases, 
namely, with either (6, 3, +i) or (8, 2, +i) in the desert, 
fully consistent solutions do exist. Moreover, in both 
these settings, the upper limit on the B — L scale is 
pushed up by several orders of magnitude with respect 
to the naive estimate based on the extended survival hy- 
pothesis [ll|, thus opening a room for a natural imple- 
mentation of the renormalizable seesaw in this class of 
5O(10) GUTs. 

1. Consistent setting 1: light (6, 3,+|) 

a. Stable vacua with a light (6,3,+^): Let us first 
quote the results of a dedicated numerical scan of the 
parametric space aiming at the identification of the stable 
vacua supporting an accidentally light (6, 3, +§)• Confin- 
ing all the dimensionless couplings into the conservative 
[— 1,+1] range stable vacua are confined in the domain 
wbl > 0, < 0, /3 4 > 0, a > -0.1 and | 72 | < 0.6. 

b. One-loop unification with a light (6, 3, +-|)-' Sam- 
ple regions of the parametric space that support a consis- 
tent scalar spectrum and, at the same time, provide a vi- 
able gauge-coupling unification, are depicted in FIGs.[TJ[2] 
andO It is worth noticing that: i) In the fully consistent 
settings (6, 3, +|) is pinned down to a relatively narrow 
region around 10 11 GeV, cf. FIG.[TJ ii) For all consistent 
configurations we find \uir\ <C ujbl so these settings gen- 
erally prefer an intermediate Z c 2l2r\bl stage, cf. FIG. [2] 
Hi) There is just a little room left if the current Super- 
K limits get improved considerably in future, see FIG. [T] 
Moreover, since the two-loop effects tend to further lower 
the GUT scale with respect to the one-loop estimates 
(even as much as half an order of magnitude) fill ] . this 



FIG. 1. M(6, 3, +i) — ujbl correlation in the case of a light 
(6, 3, +-|) multiplet in the desert. Various levels of gray cor- 
respond to domains accessible for different GUT-scale limits, 
cf. Section Till A 21 M(6, 3, +^) can vary only over a couple of 
orders of magnitude (for the current SK limit) and the range 
is likely to shrink considerably in future. 



u BL [GeV] 




FIG. 2. \ujr\ — ujbl correlation in the case of a light (6, 3, +|) 
multiplet in the desert. The color code is the same as before, 
cf. Section [HI A 21 In all of the allowed region \ujr\ <C ujbl so 
this setting prefers an intermediate 3 c 2_l2_r1sl stage. 



class of scenarios may be ultimately testable at Hyper-K. 
iv) Finally, the actual upper limit on the B — L scale is 
stretched to almost 10 15 GeV and it slowly decreases for 
stronger proton-decay limits, cf. FIG. [3] 

However, one should be more careful here because 
these results can be biased by the stability of the numer- 
ical approach we are using, cf. Section IIII B 21 Namely, 
the system of equations implementing the unification 
constraints can be efficiently solved for the position of 
(6, 3, +|) and for the overall shift of the spectrum if and 
only if (6, 3, +|) is considerably lighter than the next-to- 
lightest threshold at play (typically a gauge boson associ- 
ated to the 2f>1bl — > ly breaking); otherwise it becomes 
highly non-linear and, hence, difficult to handle. How- 
ever, as one can see in FIG.|H for the estimate of the up- 
per limit on a this issue is less important because some of 
the couplings (namely, /?4 and /3 4 ) turn non-perturbativc 
yet before this issue really occurs. Moreover, the shape 
of the new upper limit on the B — L scale is such that 
one is likely to miss solutions in the lower- B — L regime 
which is not of the utmost importance here. 
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FIG. 3. \u)r\ — \a\ correlation in the case of a light (6, 3, +|) 
multiplet in the desert. The color code is the same as before, 
cf. Section IIII A 21 B — L as high as almost 10 15 GeV can 
be reached for the current SK limit, with the best Hyper-K 
sensitivity limit the maximum lowers to fewxl0 14 GeV. 
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FIG. 4. Correlation between the mass of the (6, 3, +^) thresh- 
old and the allowed B—L scale a. There are two basic stability 
issues that bias the estimate of the span of the allowed do- 
mains: first, there is the technical requirement we impose on 
the hierarchy between the lightest and next-to-lightest thresh- 
olds, i.e, (6, 3, +|) and the gauge sector associated to the 
3 c 2_l2_r1_bl breaking which cuts the parametric space from 
below right; for large <r's this, however, becomes irrelevant 
because some of the couplings (namely, /?4 and become 
non-perturbative yet before such a numerical instability can 
affect the relevant upper bound. 



c. A specific example with a light (6,3,+^).' The 
"effective" SM gauge coupling evolution with a light 
(6, 3, +^) is exemplified in FIG. [5]where the values of the 
input parameters as specified in the left row of TABLE ITO1 
have been used and r is calculated so that the desired 
M(6,3,+|) = 5.57 x 10 11 GeV is obtained. Note that 
the small I72I region turns out to be preferred for larger 
values of \a\ and that we have chosen a solution with 
relatively small | A4 1 and A4 just to optically improve the 
expected "clustering" of the (3, 2, +|) and (3, 2, +i) mul- 
tiplcts at around 10 15 GeV (cf. FIG. [5]) due to their com- 
mon origin within (3, 2, 2, +|) of Z c 2iJIr\bl- A more de- 
tailed information about the relevant bosonic spectrum 
underlying the gauge unification in this setting is given 
in TABLE HV] of Appendix [D] 
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FIG. 5. Unification of the effective SM gauge couplings in 
a sample setting with a light (6,3, +|) multiplet (here at 
around 5.6 x 10 11 GeV, cf. Section [Til C 1 cfl with the shaded 
area magnified on the lower panel. The short 3 c 2_l2_r1sl 
stage is clearly visible here. The small circles indicate the 
positions of various thresholds (for details, see TABLE IIV[I 
inflicting changes in the three curve's slopes. The almost- 
vertical solid and dashed lines correspond to the current and 
future proton-decay limits, cf. Section [III A 21 The displayed 
setting is compatible (at one loop) with the current SK limit, 
but it can be refuted by the Hyper-Kamiokande. The dotted 
vertical line indicates the position of the B — L scale. 



2. Consistent setting 2: light (8, 2,+i) 

a. Stable vacua with a light (8, 2,+i): Turning our 
attention to the remaining option of a light (8,2,+^) 
in the desert it is possible to show that (for all dimcn- 
sionlcss couplings smaller than 1 in absolute value) there 
are always tachyons in the scalar spectrum outside the 
following domain: \(jJbl\ < l^flj, P'± < 0, ao > —0.05, 
I72I < 0.6, I72I < — 0.8(3' 4 . Moreover, only one of the 
eigenstates of the mass matrix (|B11[) can be consistently 
lowered. 

b. One-loop unification with a light (8, 2, +\): Sam- 
ple regions of the parametric space that support a con- 
sistent scalar spectrum and, at the same time, provide a 
viable gauge coupling unification, arc depicted in FIGs.|6l 
[7]and[8l Note in particular that: i) In the fully con- 
sistent settings (8,2, +|) is narrowed down to the lower 
part of the desert (it is always below 10 10 GeV) and even 
more so if proton lifetime limits get improved in the near 
future, cf. FIG. |H1 Nevertheless, the allowed domain for 
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TABLE III. Parameters underpinning the two sample settings 
with a light (6, 3, + §) (left) and a light (8, 2, +±) (right), re- 
spectively. The value of the r parameter can be obtained from 
the requirement that the relevant light threshold has the mass 
specified in the last row. 
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FIG. 7. \lobl I — |wr| correlation in the case of a light (8, 2, + |) 
multiplet in the desert. Various levels of gray correspond 
to domains accessible for different GUT-scale limits, cf. Sec- 
tion [TTT^2] In the whole allowed region \uj bl\ <C 0Jr so this 
setting always exhibits an intermediate 4c2l1r stage. 
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FIG. 8. An interesting \ljbl\ — correlation in the case 
of a light (8,2,+^) multiplet in the desert. The color code 
is the same as before, cf. Section IIII A 21 B — L as high as 
10 14 GeV can be reached and, remarkably enough, unlike in 
the (6,3,+^) case the maximum reach is insensitive to the 
proton lifetime limit. 



FIG. 6. M(8, 2, +|) — \ur\ correlation in the case of a light 
(8, 2, + |) multiplet in the desert. The color code is the same 
as before, cf. Section fill A 21 M(8, 2, +5) can vary over many 
orders of magnitude in the lower part of the desert, and it is 
pushed down for increasing proton lifetime. 



the light threshold is much wider than in the previous 
case, see FIG. Q] hence, this scenario is likely to be more 
robust to the changes inflicted by two-loop effects. In the 
extreme case this class of models requires (8, 2, + i) close 
to the EW scale with possibly interesting collider effects. 
ii) For all consistent configurations we find \u)bl \ <C \cjr\ 
so these settings generally prefer an intermediate 4c2l1r 
stage, cf. FIG. [7] Hi) The upper limit on the B — L scale 
is pushed to about 10 14 GeV and, unlike in the previous 
case, it is rather insensitive to possible future improve- 
ments of the proton lifetime limits, cf. FIG. [5] 

Note also that there is no problem with the numer- 
ical stability here, cf. Section IIII B 21 because the gap 
between the mass of (8, 2, +5) and the next-to- lightest 
threshold preferred by the unification constraints is huge. 



c. A specific example with a light (8,2,+i): The 
"effective" SM gauge coupling evolution with a light 
(8,2,+i) is exemplified in FIG. |H] where the values of 
the input parameters as specified in the right column 
of TABLE Mil have been used and r is fixed so that 
M(8,2,+±) = 2.3 x 10 4 GeV. Note that the small | 72 | 
region turns out to be preferred for larger values of |c| 
and we have chosen smaller | A4 1 and A4 just to opti- 
cally improve the expected "clustering" of the remnant 
of the (15,2, +|) multiplet oi A^c^lAr (where the light 
(8, 2, +±) comes from) at around 3 x 10 14 GeV, cf. FIG.|U 
A more detailed information about the relevant bosonic 
spectrum underlying the gauge unification in this setting 
is given in TABLE [V] of Appendix |DJ 



3. Further remarks 

It is perhaps interesting to note that neither (6, 3, +3) 
nor (8, 2, +5), even if pushed well below the GUT scale, 
are any way problematic for cither BBN or d > 6 proton- 
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FIG. 9. The same as in FIG. here with a light (8,2,+|) 
multiplet (at around 2.3 x 10 4 GeV, cf. Section IllTTjTcl) in- 
stead of (6, 3, + §); for details of the spectrum see TABLE [VI 
The short 4c2l1r stage as well as the "fake unification" fea- 
ture is clearly visible here. The displayed setting is compatible 
(at one-loop) with the current SK as well as possible future 
Hyper-Kamiokande limits. 



decay, cf. Sects. IIII A 2 bl and MI A 31 As for the former, 
their direct coupling to quarks and leptons (by means 
of the 16^16^126^ SO(10) Yukawa interaction) makes 
both of them decay fast enough via a simple tree-level 
diagram. Similarly, neither (6, 3, +^) nor (8,2,+i) can 
mediate the effective d = 7 (B + L conserving) proton 
decay. 



IV. THE MINIMAL SO(W) GUT REVIVED 

The previous analysis reveals several regions of the 
parametric space of the non-SUSY 50(10) Higgs model 
based on the reducible representation 45 h © 126// that 
can consistently support 5*0(10) — > SM symmetry break- 
ing chains compatible with the electroweak data and 
the current proton decay limits and, simultaneously, 
admit for a large-enough B — L breaking scale for 
a natural implementation of a rcnormalizable seesaw. 
Hence, this simple Higgs model is ready to be upgraded 
to a full-featured, potentially realistic and predictive 
50(10) GUT. 

In doing so, the central question to be addressed be- 
fore approaching any of the ultimate goals of such a pro- 
gramme (e.g., a detailed prediction of the proton lifetime 



and the relevant branching ratios) is the structure of the 
Yukawa sector. 



A. Yukawa sector of the minimal SO(W) GUTs 

It is easy to see that the Higgs model containing just 
45// and 126// can not, at rcnormalizable level, support 
a viable Yukawa sector as there is only one contrac- 
tion available in such a case, namely, 16f/ 126 16_f126^. 
Hence, the flavour structure is entirely governed by a sin- 
gle (symmetric) matrix of Yukawa couplings / 126 and no 
mixing nor featured fermionic spectra can be generated. 

The minimal potentially realistic extension of the 
45// © 126// setting amounts to adding an extra 10- or 
120-dimensional representation which can smear the de- 
generacy of the effective Yukawa matrices across differ- 
ent fermionic species; for a more detailed discussion see, 
e.g., [HI or, more recently, |l6|. In this respect, it is 
interesting to quote namely the results of the new nu- 
merical analysis [HI attempting to fit the SM flavour 
structure onto the effective mass matrices emerging in 
both the 126// © 10// as well as the 126// © 120// cases: 
Interestingly, the former option is strongly preferred and, 
moreover, successful fits require a dominance of the type- 
I seesaw contribution^. However, as interesting as these 
results are, they are still not entirely decisive as there 
are various sources of uncertainties^ that have not been 
taken into account in [48| . 

Nevertheless, the Higgs sector based on 45// ffi 10// © 
126// is clearly the first choice; not only it has a better 
chance to be compatible with the fermionic data, but the 
addition of an extra 10// rather than a larger multiplet 
like 120// only minimally disturbs the results obtained in 
the previous parts, see also the comments in Sect. IIII B 21 

For the sake of completeness, let us reiterate the 
Yukawa-sector sum-rules relevant to this setting. In full 
generality, one can write a renormalizable Lagrangian 
densitjO 



C 3 16 F (/ lu 10// + / 2 lu 10^ + / i2b 126^)16 F + h.c, (17) 

which is parametrized by three complex symmetric ma- 
trices fl° 2 and / 126 . It leads to the following general 
tree-level formulae for the effective SM quark and lepton 



1,3 This feature is closely related to the need to avoid the b-r Yukawa 
unification in the non-SUSY settings which, however, is generi- 
cally favoured by type-II seesaw. 

14 In particular: i) the weights of the SM-doublet VEVs entering 
the relevant sum-rules, cf. Eqs. 1181 1. were taken uncorrelated, 
ii) the running fermionic masses were extrapolated to the GUT- 
scale vicinity under the bold assumption of no thresholds in the 
desert and Hi) higher order corrections to the relevant sum-rules 
inherent to non-SUSY settings were not taken into account. 

15 Note that 10* of 5O(10) is equivalent (in the representation 
sense) to 10, so both 16 2 10 and 16 2 10* are allowed in non-SUSY 
scenarios. 
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mass matrices^ 



M u 


= ^ V° 4 


-y 10 ^* + 


^126 126 


M d 


= Yl%f 4 


-F 2 10 ^ * + 


yl26 126 


M e 


= W 4 


-Y 2 10 vT- 


3F 126 V 126 ; 




= tfV 4 




3Y 126 ^ 26 , 


M M,I 




M M,II = 


c II y 126 w 126 



(18) 



where and Y 126 are proportional to the correspond- 
ing /-matrices in (|T7|) . the subscripts D and M denote 
the Dirac and Majorana segments of the neutrino mass 
matrix and c I,n (for type-I and type-II seesaw, respec- 
tively) combine various extra numerical factors such as 
the relevant Clebsch-Gordan coefficients. Let us note 
that, unlike in SUSY, here the 10// of SO(10) can be 
populated by real components [f|, which would further 
reduce the number of independent couplings - indeed, 
in such a case, the second term in ([!?]) would be just a 
repetition of the first one. However, this setting is patho- 
logical as it leads to a GUT-scale near-equality of the b- 
and i-quark masses, cf. [l6j . 

It should be stressed that in the full model, the 
projections of the SM Higgs VEVs onto the indi- 
cated components of the relevant SO(10) multiplets 



v? = ((l,2,+±)io), vf = ((l,2,-i) 10 ), = 
((l,2,+i)i26-) and v 1 / 6 = ((l,2,-i)i 26 .) as well as the 
VEVofthe type-II SU(2) L triplet w 126 = ((1, 3, +1) 126 «) 
are calculable functions of the parameters entering the 
scalar potential. Note also that, in full generality, the 
formulae (fT8|) correspond to the two-Higgs-doublet real- 
ization of the SM Higgs sector; assuming that only one of 
the doublets survives down to the electroweak scale (i.e., 
implementing a single fine-tuning in the relevant gener- 
alization of the mass matrix (|B9[) ) one should further 



,126 



or vf 



,126 







assume eithei0 v^ 

In connection to this, one should mention a couple of 
other interesting features inherent to the models with 
45// that have no counterpart in many other settings in- 
cluding the popular variant with the GUT symmetry bro- 
ken by 54// [3(1 HH, H9|: i) First, in the former case, a 
significant admixture of both the 10// and 126// compo- 
nents within the SM Higgs doublet, which turns out to 
be essential for realistic fits of the Yukawa system (TTS)) 
to the quark and lepton data, is naturally obtained via 
the direct mixing term 10//126^45//45//. In the latter 
case, however, there is no such a term and the mixing is 
governed namely by the 10//126^126//126// contraction 
which, however, yields and extra suppression of the or- 
der of M] 3 _ l /Mq. ii) Second, the settings in which the 



GUT symmetry is broken by the VEVs of 45// generally 
feature an almost-automatic suppression of the type-II 
seesaw which, as mentioned previously, is not only wel- 
come due to the generic GUT-scale non-equality of the b 
and r Yukawas in non-SUSY settings [5(j, but it seems to 
be even crucial for successful Yukawa fits, cf. 48j. Indeed, 
on general grounds, one expects that in theories in which 
the D-paritj{3 is broken before the SU(2)r symmetry, 
the type-II contribution to the light neutrino masses is 
naturally suppressed by a factor of M%_ l /Mq with re- 
spect to the type-I term 55j. Again, this is not the case 
in models based on the 54// where the D-parity is pre- 
served down to the S'J7(2)/j-brcaking scale and, thus, no 
extra suppression of type-II seesaw occurs. 



B. Predictivity and testability 

Concerning the predictivity of the renormalizable 
model based on 45// © 10// © 126// in the Higgs sector, 
there are several aspects worth a comment here. 

i) Yukawa sector complexity: There are in general 
three independent complex symmetric matrices enter- 
ing the effective sum-rules (|T5j). to be compared to just 
two such structures encountered in, e.g., the minimal po- 
tentially realistic Yukawa sector in supcrsymmetry [17| . 
This, however, does not necessarily imply a loss of pre- 
dictivity in the Yukawa sector: First, the weights of the 
SM VEVs entering Eqs. (TT8"| are in general stronger cor- 
related here than in the minimal SUSY case (cf. (20j ) 
because here the doublet mass matrix (|B9[) is lower- 
dimensional. Second, with only one doublet pushed down 
to the electroweak scale, the system (TTS|) is simplified and 
the correlations among different species become much 
tighter. This is also well reflected by the preliminary 
results of a dedicated numerical analysis (56|. 

ii) Vacuum stability: Unlike in the (global) SUSY case 
where the positivity of the scalar mass-squares is auto- 
matic in any SUSY-preserving vacuum, the consistency 
requirements here narrow the potentially interesting do- 
mains of the parametric space down to just few small 
patches (for instance those identified in Sect. IIII[) . On 
the other hand, given the higher number of contractions 
available in the non-SUSY case even at the renormal- 
izable level (to be compared to just several such terms 
entering the Higgs superpotential in SUSY) the set of 
SM-like vacua is clearly higher-dimensional (see the num- 
ber of parameters in TABLE IlIII versus a single complex 
parameter in SUSY, cf. 0). 



Let us remind the reader that a good grip on the Yukawa cou- 
plings is also necessary for a reliable account of the d = 6 proton- 
decay amplitudes because they depend on the matrix elements 
of the unitary transformation bringing the quarks and leptons 
from the current to the mass basis. 

Actually, only the former option has a chance to work in prac- 
tice because the latter immediately implies an apparently wrong 
relation M d = M( so we discard it. 



D-parity is a discrete symmetry actin g as a charge conjugation 
in the left-right symmetric context [51],[52| , and, as such, it plays 
the role of a left-right symmetry (enforcing, for instance, equal 
SU(2) L and SU(2) R gauge couplings). As a part of the SO(10) 
algebra, it is a good symmetry until it gets broken by either a 15- 
odd singlet in 45// (or 210//) or by any S(7(2) fi -breaking_yEV, 
thus allowing for a left-right asymmetric scalar spectrum l53tl54l| . 
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Hi) Radiative corrections: Unlike many popular SUSY 
50(10) variants, cf. 0, H3] the model under considera- 
tion is asymptotically free (with b = —12) so it remains 
weakly coupled up to the Planck scale. On the other 
hand, its radiative structure is much more involved than 
that of the simplest SUSY scenarios and consistent cal- 
culations are technically much more demanding. In this 
respect let us reiterate that the nature of the problem 
calls for a two-loop RGE analysis based on a detailed 
knowledge of the one-loop spectrum and in this work we 
have just performed the first steps in that direction. 

Hence, without a detailed analysis it is rather difficult 
to assess the predictive power of the model under consid- 
eration. Nevertheless, even the first results obtained in 
Sect. lIIII indicate that the up-coming large volume experi- 
ments such as Hyper-K can impose very strong cuts to its 
(already rather constrained) parametric space, possibly 
covering the entire remaining volume. 

Sometimes, it is suggested to further enhance the 
Yukawa-sector predictivity of the non-SUSY models by 
imposing an extra global U(l) symmetry of the Peccei- 
Quinn (PQ) type [58|, [59[ which, if it transforms 10# 
non-trivially, forbids one of the f 10 2 couplings in the 
Lagrangian (fT7|) . Since, in that case, also 126h would 
have to be PQ-charged, such a symmetry would be bro- 
ken at the same scale as U(1)b-l, thus linking the PQ 
symmetry-breaking scale to the neutrino masses. In this 
respect, it is very interesting that a seesaw scale in the 
preferred 10 13 ^ 14 GeV ballpark is indeed very close to the 
10 9H ~ 12 GeV PQ-symmetry-breaking window favoured by 
astrophysics and cosmology (see e.g. and there are 
several attempts in the literature to construct a viable 
unified model along these lines (see, e.g, [6l| . or more 
recently [l6|). On the other hand, since (126 h) can not 
break the rank of 5*0(10) <g> U(1)pq by more than a sin- 
gle unit, a global linear combination of U(1)pq, U(1)r 
and U(1)b-l survives down to the electroweak scale and 
only there it gets finally broken by the electroweak dou- 
blet^); this, however, is unacceptable as the EW-scale 
PQ-symmetry breaking gives rise to an easily visible ax- 
ion (62l. [63j. Thus, a consistent implementation of this 
interesting scheme calls for a further complication of the 
Higgs sector, which we shall not entertain here. 



V. CONCLUSIONS AND OUTLOOK 

In this work we have been concerned with a class of 
simple renormalizable S*O(10) Higgs models in which the 
first stage of symmetry breaking is triggered by the 45- 
dimensional adjoint Higgs representation. These set- 
tings, discarded a long ago due to inherent tree-level 
tachyonic instabilities developing in most of the physi- 
cally interesting scenarios, have been recently shown to 
be revived by quantum effects. However, many impor- 
tant aspects of these scenarios, such as, for instance, the 
allowed ranges for the unification as well as various inter- 



mediate scales, as important as these are for any realistic 
model building, were never studied in sufficient detail 
to allow for a qualified assessment of their physical rele- 
vance. 

Focusing on the variant with 45# © 126# in the Higgs 
sector, we worked out the complete tree-level spectrum 
and, with such an extra information at hand, performed 
a simple analysis of the gauge-unification patterns. Un- 
like in the previous studies based on the minimal-survival 
hypothesis, that show a no- go for scenarios with the 
B — L scale above 10 10 GeV [11(, we found several do- 
mains in the parametric space of these models that can 
support a consistent gauge-unification with B — L as 
high as 10 14 GeV without encountering any tachyonic 
instabilities or proton lifetime issues. The key to this 
unexpected behaviour is an accidentally light threshold 
in the desert which affects the gauge-unification picture 
in a suitable way. We identified two distinct classes of 
such viable solutions: in the first case, an intermediate- 
scale multiplet transforming as (6, 3, +-|) of the SM sup- 
ports S*O(10) — > SM descents featuring a short St/(3) C (E> 
SU(2) L <g> SU{2) R ® U(1) B -l intermediate-symmetry 
stage, while the second option including a relatively 
light (8, 2,+i) supports S*O(10) breaking chains pass- 
ing through the SU(A) C <g> SU{2) L (g) U(1) R intermediate 
symmetry. Remarkably enough, in all the cases of in- 
terest, the unification scale turns out to be rather close 
to the current Super-Kamiokande proton-lifetime lower 
bound. 

This, however, opens up an intriguing possibility to 
construct a simple, renormalizable and testable 15*0(10) 
GUT with 45H©126_f/©10 ff in the Higgs sector which, in 
view of the recent failure of the simplest supcrsymmetric 
5*0(10) model [26|, [27|, can even be viewed as the new 
minimal potentially realistic S*O(10) GUT. 

Nevertheless, this study provides only the first glimpse 
on the ultimate viability of such a framework and there 
is much more still to be done. Let us reiterate that sim- 
ple non-SUSY models suffering from significant tree-level 
vacuum instabilities generically call for a refined two-loop 
RGE approach (assuming one-loop scalar spectrum) be- 
cause only in such a case the tachyons are really under 
control. In this respect, the results of the current analy- 
sis, taking into account only the minimal set of radiative 
corrections necessary for the scalar spectrum regulariza- 
tion, can quantitatively (though not qualitatively) differ 
from those to be obtained in a future full one-loop effec- 
tive potential analysis. 

Remarkably enough, extrapolating the relative size 
and direction of the two-loop effects observed in [ll| to 
the current scheme, the chances for its ultimate testa- 
bility at future experiments look rather promising. In- 
deed, any further significant decrease of the maximum 
allowed unification scale due to two-loop effects would al- 
low the up-coming large- volume facilities such as Hyper- 
Kamiokande to scan over the full physically interesting 
domain in the parametric space of this class of models. 
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</>) 2 (EE*) 2 = (j>ijHklmni^'klmnj 
#)o(SS*)o = <j>ij<j>ij^klmno^*kl m no 
<W>) 4 (EE*) 4 = (j)ij4>kl^mnoij^m n okl 



4<(E£*)4' = (l)ij(j)klT, mno ik'S^ nno ji 



i(EE) 2 



&ij (Pik ^Imnoj ^Imnok 



(A7) 



2 (E*E*) 2 = 0ii0ifeS* TOnoj S* mnoi . . 

We have checked that this constitutes a complete set of 
50(10) invariants for the 45-126 system at the renormal- 
izable level. 



Appendix A: Details of the scalar potential 

The 2-index and 5-index completely antisymmetric 
tensors of 50(10) are labelled respectively as cf>ij and 
4>ijkim ■ Given the dual map 



(Pijklm = ~-^ e ijklm,nopqr4 ) n,opqr i 

(Al) 



we can define the self-dual and the anti-self-dual irre- 
ducible components of 4>ijkim as 

^ijklm = ^= (j>ijklm + ^ijklni^j , (A2) 



^ijklm — (fiijklm 4>ijklm 



(A3) 



Then the relevant contractions in the scalar potential 
of Eq. ([1]) are given by 



)0 = <pij<pij, (SS*)o = ^ijklm^ijklm ( A4 ) 
)o(<W>)o = (/>ij<f>ij<f>M<f>kt (A5) 
) 2 (</>0) 2 = <f>ij(t>ik(i>lj(t>lk 

(EE*) (EE*) = ^'ijklm'^ijklm^'nopqr'^'nopqr 
(EE*) 2 (EE*) 2 = Y.ijklm^ljkln^opqrmX'opqrn 
(EE*)4(SE*)4 = Tjijklm^ljkno^pqrlrnJ^pqrno (A6) 

(SE*) 4 /(EE*) 4 ' = Y; ij ki m T<* jkno T, pqr i n Y,* qrmo 

(EE) 2 (EE) 2 = ^ijklrn^ijkln^-'opqrrn^-'opqrri 

/v^*x^*^ /V , *V*\ y~i* v^* V* V* 

[Zj Zj )2{2-< Zj )2 — 2^ i j klm n i: j kln n opqrm l Jopqrn 



Appendix B: The scalar spectrum (45.tr © 126h) 



1. Vacuum manifold and stationarity conditions 

The scalar potential Eq. (JXJ) evaluated on the SM vac- 
uum parametrized by ubl, &> r and a, cf. Eq. ([5]) reads 

(V)= -m 2 (3w| L + 2w|) + a (12lj 2 r uj% l + 9u>% L + 4u%) 



•f{3u% L + 2u>%)-2v*\*\'+4X \a\ 



-2t (3u) B l + 2u R ) \a\ 2 + 2a (3w 



2 

BL 



2col) \a\ 2 



-4/?4 (6uj r oj bl + 3oj bl + uj 2 r ) \a\ 



(Bl) 



It is perhaps worth noting that not all couplings in ex- 
pressions ©-(HI) are present here; the reason is the ab- 
sence of suitable terms quartic in the available VEVs in 
some of the contractions. As an example consider 7J2, 
the coefficient of the 126 B contraction in Eq. ([3]), which 
enters neither the vacuum manifold nor the stationary 
conditions or the tree-level spectrum. This can be under- 
stood by looking at the decomposition of the the relevant 
invariant under SU(5)<g>U (1) z which never contains more 
than a single submultiplct (1, +10) that is the only com- 
ponent of 126 h that can receive a SM-preserving VEV 
(recall that (126^) preserves SU(5)). Indeed, one has at 
best 126^ D (l,+10)(T5,-6)(T5,-6)(50,+2), i.e., three 
derivatives are needed in order for r\2 to enter anywhere. 
A similar reasoning can be applied to the other couplings, 
hence fully justifying their presence/absence within all 
the relevant structures. However, most of such couplings 
(e.g., A 2 , A4 etc.) reappear in the tree-level broken-phase 
mass matrices and, ultimately, all of them appear at the 
full one-loop effective potential level. 

The corresponding stationary equations can be conve- 
niently rewritten as 
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d(V) 3d (V) 



-/i 2 + a (6w|x, + 4u%) + a 2 (urUbl + w% L + u%) + 2a\<r\ 2 + 2/3' 4 \a\ 2 ,(B2) 



4(ujbl - UJr) \3 duiBL 



2 d(V) 
Mr-. 



ubl- 



d(V) 
da 



2a* 



Oujr 



a 2 u} R oj B L ( w bl + w B ) - r\a\ 2 + 2j3' 4 (3uj B l + 2oj r ) \a\ 2 , 



4X \a\ 2 + r (3u BL + 2uj r ) + a (3uj 2 bl + 2ui%) - 2(3' 4 (6uj r uj bl + 3u 2 BL 



(B3) 
(B4) 



r 



which hold away from the standard SU(5) ® U(l)z vac- 
uum (UJBL = Wr). 



in this sector one can identify 12 Goldstone boson modes 
with 



2. Tree-level scalar spectrum in the SM limit 



M 2 (3,2,- 



r2/Q o 5>\ q 



(B7) 



Let us label the scalar states with respect to the 
SU(3) C <g) SU(2) L ® t/(l)y algebra. Applying first the 
stationary conditions in Eq. (|B2[) and Eq. (|B4|) one finds: 



The remaining components of 45 h rnix with those in 
126# and will be discussed below. 



a. States with components in 45h only 



Besides the classical pair of the would-be tachyons with 
mass formulae 

M 2 (l, 3, 0) = -2a 2 {uo BL - cJr)(uj bl + 2cj r ) , (B5) 
M 2 (8, 1, 0) = -2a 2 (u) R - u BL ){ur + 2lo bl ) , (B6) 



b. States with components in 12&h only 



Starting with pure states with components in 126 B one 
has: 



M 2 (l,l,+2): 

M 2 (l,3,-1): 
M 2 (3,l,+§) 
M 2 (3,3,-i) 
M 2 (6,3,+i) 

M 2 (B,l,-f) 
M 2 (6,l,-i) 
M 2 (6,l,+§) 



3A 4 
3A 4 
3A 4 
3A 4 



16A 4 )|<7| 2 -4^(t-6/3^ sl ) , 
2\' A ) \a\ 2 - 2 {3u BL + w R ) (r - 2^ R ) 



2(4A 2 
8(2A 2 

= 4 (3A 2 + 3A 4 + 4\' A ) \a\ 2 + 2 (oj B l {0 a (ubl + 2w R ) + P a w B l) - r (uj B l + 2u R )) , 
= 4 (3A 2 



4A^) \a\ 2 - 2t (2u BL + u R ) + 2/3 4 w| L + 4(3' 4 {3uj r uj bl + 2lu 



BL 



2 (4A 2 + 3A 4 + 16Ai) W? + 2 (^bl + u R ) (2ft {2u BL + w R ) - r) , 
2 (4A 2 + 3A 4 + 16Ai) \a\ 2 + 4u BL (2# (^bl + 2u R ) - t) , 



4(3A 2 
8(2A 2 



3A 4 
3A 4 



■ 4A 4 ) \a\ 2 ~ 2t (2uj bl + u R ) + 2/3 4 lo 2 r + 4(3' A (3u R u; BL + 2lo 2 bl + oj 2 r ) 
• 2X' A ) \a\ 2 - 4 (ujbl + u R ) (t - 2#w Bi ) , 



(B8) 



r 



while those developing higher-dimensional mass matrices 
are: 



M 2 (l,2, 



4(3A 2 



3A 4 + 4A^)|ct| 2 - 3r {uj BL + w R ) + {-4uj r u bl + 7uj 2 bl + u%) + 3^ 4 (4lo r lu bl + 3lu 2 bl + u%) 



27l 



12 \ W R - ^BL 

272 {u% - w% L ) 
2(4A 2 + 3A 4 - 8A^)|a| 2 - r (3uj bl + u R ) + f - (4lu r lo B l + 7uj 



2 

BL 



3^ 4 (4u> R L> BL 



(B9) 
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w2/„ o , 7n _ ( 2 ( 4A 2 + 3A 4 + 16A 4 ) \v\ 2 - r (lu bl + 3lj r ) + {lo bl - u R ) 2 + fa 4 (16oj r lj bl + 5u% L + 3u%) 



2 72 (w| 



BZ, 



272 (w| - w^) 

4 (3A 2 + 3A 4 + 4A^) \a\ 2 - t {5uj bl + ur) + ^ (lu bl + w R ) 2 + fa (l6u R co BL + hu 2 BL + 3tu 2 R ) 



(BIO) 



M 2 (8,2, 



1 / 4 (3A 2 + 3A 4 + 4A^) M 2 + ± + (6$ (3w B l + w«) + fa (u BL + w«) - 6r) 
~ 1] \ 2 72 *(a,|-c4 L ) 

272 (w| - w| L ) 

2 (4A 2 + 3A 4 + 16A 4 ) |<r| 2 - r (3w BL + ojr) + & (lu bl - lor) 2 + 3/3 4 (4u r oj bl + 3lu 2 bl + u 2 R ) 



(Bll) 



2(4A 2 + 3A 4 + 8A 4 )|ct| 2 - 2r (w B£ + lj r ) + fa (u% L + uj 2 r ) + Afa 4 (3uj r uj bl + 2u 2 BL + u)%) 
M 2 (3,l,+±)=| 4rt(u% L -<4) 

2 % /2(-8A 4 |a| 2 +P4.WRUB1) 
472 (w|i - w 2 ,) 

4 (3A 2 + 3A 4 + 4A 4 ) |a| 2 - 2r (2w BL + w fl ) + fa {uj 2 bl + w 2 R ) + 4# (3w fl w BL + 2w| £ + w|) 



2 x /2(-8A 4 |a| 2 + ( 4 ^ BL ) \ 

, (B12) 

2(4A 2 + 3A 4 )|a| 2 + 2 (2 fa (3w fl Wfl £ + 2u% L + lo 2 r ) + fa (lo 2 bl +u r )-t (lu bl + lur)) J 



where the mass matrices above are spanned on the follow- 
ing bases (listing only the column basis vectors; the rows 
are just their conjugates): {(1, 2, -±) E », (1, 2, -|) s }, 
{(3,2,-£) s ,,(3,2,-|) E }, {( 8 ,2,-i) E ,,(8,2,-i) s } 
and {(3, l, -|) s ,, (3, 1,-1)^,(3,1,-1)1,}. For poten- 
tially ambiguous cases we use superscripts to indicate 



r 



the SU (2)r origin of the relevant components. 



c. Mixed states with components in both 45h and 126.tr 

Finally, the remaining components of 126# that mix 
with those in 45/f are listed below. 



Af 2 (l 1 +1) = I 2 (~ a ^ BL ^ BL + ^ + ^ - 2 ^) M 2 ) 2 ( 2 / ? 4 (3wsl + Wfl) + fauR -r)a*\ 
1 2 ( 2 # (^bl + u>r) + fauR -t)o 2lo r (2 fa" {3lo bl + cjr) + faujR - r) ) ' 



M 2 f3 1 [ 2 (-«2Wii(w B i + + G84 - 2 /? 4 ) kl 2 ) " 2 (4/? 4 (w BL + w H ) + fauj BL -r)a\ 

3 V -2 (A fa (u BL + ur) + fau BL - t) a* 2u BL (A fa {lo bl + u R ) + fau BL - r) j ' 

M 2 (3,2,+I) = 

/ 2 (-2a 2 u BL u R + {fa - 2fa 4 ) \a\ 2 ) (2fa 4 (hw BL + 3oj R ) + fa (u B l + ur) ~ 2r) a* 

(2 fa (5uj bl + 3u> R ) + fa (lu bl + lur) -2t)<j 2 (lu bl + lur) (2fa 4 (5lu bl + 3lur) + fa (uj bl + w R ) - 2r) 

\ 4 7 ^ (UJR - UJ BL )<7 2 72 * (UJ 2 R -UJ 2 BL ) 

472 (ur - lo B l) ct* 
27 2 K-a;| L ) | , (B15) 

8 (2A 2 + 3A 4 + 2A 4 ) \<j\ 2 - r (5wij L + Slur) + & (lu bl -lu r ) 2 + fa 4 (8urlu bl + 5lu 2 bl + 3u%) 
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M 2 (1,1,0) 



( 2 (l2a cj| L + a 2 (uJ B L - u R )(2u BL + w R ) - 6/? 4 |cr 

4V6(2a wsLWi?-^kl 2 ) 
v/6 (t + 2aw BL - 4/34 (w B£ + w K )) ff* 
V6 (r + lauiBL - 4/3 4 (w BL + w fl )) cr 

V6(r + 2aw BL -4/?4(a; BL 
2 (r + 2aw fl - 2/3 4 (3uj BL 4 
4Aok| 2 
4A cr 2 



4^/6 (2a oWi3L ^-^|a| 2 ) 
2 (8a w| - a 2 (wBL - uj r )(uj B l + 2u R ) - 4^|cr| 2 ) 



2(r- 
2(t 

hWfi))(T \/6(t 
wr))ct 2(t4 



■ 2aw fl - 2/34 (3w BL + w R )) cr* 
- 2aw_R - 2/34 (3wbl + Wfl)) cr 
I- 2auj BL - 4/3 4 (lobl +oj r ))<j* ^ 
2aw fl - 2/3 4 (3wbl + cr* 

4A cr* 2 

4A |cr| 2 / 



(B16) 



Here the relevant bases are: {(3, 1, — |)^, (3, 1, — } 
and {(3,2,-I) ,(3,2,-I) s .,(3,2,-A) s } for the 
coloured triplets and {(1, 1, — 1)^, (1, 1, — l)s*} and 

{(1,1,0)^,(1,1,0)1,(1,1,0)^,(1,1,0)2} for the sin- 

glets. Wherever ambiguous, the superscripts denote the 
SU(4)c origin of the relevant components. 

Finally, implementing the remaining stationary condi- 
tion in Eq. (|B3| (substituting for 02) one obtains 

RankM 2 (l, 1,+1) = 1, 
RankM 2 (3,l,+§) = 1, (B17) 
RankAf 2 (3,2,+±) = 2, 
Rank Af 2 (l, 1,0) = 3. 

Together with Af 2 (3, 2, -|) = 0, cf. Eq. (|B"7) . we account 
for exactly 33 WBG bosons corresponding to the coset 
SO(10)/SU(3) C <g> SU(2) L 8> U(l) Y . 



3. Basic consistency checks 

In order to crosscheck the results given above, we shall 
now study the scalar spectrum of Appendix IB 21 in three 
physically interesting limits. In each case one should ob- 
serve a proper re-clustering of the SM multiplets accord- 
ing to the enhanced symmetry as well as extra WGBs. 

a. The flipped 5' lz> limit, ur = —ujbl 7^ and a = 

Labelling the scalar states according to the flipped 
5' lz' algebra, cf. TABLE HI the 45# components cluster 
as follows 

M 2 (24,0) = 4a 2 u; 2 , 

M 2 (10, -4) = 0, (B18) 
M 2 (1,0) = 4(10a + a 2 )w 2 . 



Similarly, for the 126# components we get: 

M 2 (l, +10) = -v 2 + 5 ((a - 4/3 4 ) u-r)u, 
M 2 (5,+2) = -^ 2 + ((5a + 6/3 4 + 4/3 4 )w-T)w, 
M 2 (10, +6) = -v 2 + ((5a + 2/3 4 - 4/3 4 ) u-3t)u, 
M 2 (15, -6) = -v 2 + ((5a - 4/3 4 ) u + 3r) ui , (B19) 
M 2 (45, -2) = v 2 + ((5a + 2/3 4 + 4/3 4 ) u + r) u , 
Af 2 (50, +2) = -v 2 + ((5a + 4/3 4 ) to - r) to . 

As expected, there are 45 - 25 = 20 WGBs. 



b. The 3 c 2l2r1bl limit, lobl 7^ 0, ujr = and a = 

Labelling the scalar states according to the 3 c 2l2 r \ bl 
algebra, the A5 B components cluster as follows: 

M 2 (l,3,l,0) = -2a 2 c4 L , 

M 2 (l,l,3,0) = -2a 2 c4 L , 

M 2 (8,l,l,0)=4a 2W | L , (B20) 

M 2 (3,2,2,-i) -0, 

M 2 (3,l,l,-|)=0, 

M 2 (l, 1,1,0) =4(6oo + <to)u% L . 

Analogously, for the 126jj components we get 

M 2 (l, 3, 1, -1) = -v 2 + 3 ((a - 2f3' 4 ) uj B l - r) uj bl , 
M 2 (l, 1, 3, +1) = -v 2 + 3 ((a - 2f3' 4 ) lu bl + r) lu bl , 
M 2 (3,3, 1,-|) = -r/ 2 

+ ((3a + 2 (/3 4 + #)) cj sl - r) w BL , 
M 2 (3,l,3,+i) = -v 2 (B21) 

+ ((3a + 2 (/3 4 + /3 4 )) w BL + t) lu bl , 
M 2 (6, 3, 1, +±) = -i^ 2 + ((3a + 2/3 4 ) w BL + r) uj bl , 
M 2 (6, 1, 3, -|) = -v 2 + ((3a + 2/? 4 ) uj bl - r) w B l , 



18 



M*(l,2,2,0) = f-^ + i^ + ^ + eflJc-ix 2 1 m2 ), (B22) 

M 2 (8,2,2 I 0)=f^ 2 + ^ 6a + /34 2 +6 ^^ 2 w - 2 ^- M2 V (B23) 

«/2/o i i _ f -^ 2 + ((3a + /34 + 2^)w B L + r)w BL 4 72 w 2 L \ m „.x 

( ' ' ' +3 J " ^ 4 72 * W | L + ((3a + ^ + 2fi) u BL r) u BL ) ' ^ 

J ^ ( 3 >2) 2, + |)=f- 1/a + ' ((6a + ^- 2 f* )WM +4r) ^ 2 t „ - 272W - X , V(B25) 



where the matrices above are spanned over c. The Ac2l1r limit, ur 0, ujbl = and a = 
{(3,1,1 -i) E . ) (3,l,l-^)s},{(3 ) 2,2-f) E .,(3,2 ) 2-f) E } ) 

{(l,2,2,0) s », (1,2,2,0) E } and {(8,2,2,0) E » , (8,2,2,0) s }. Again, as anticipated, the clustering of the scalar spec- 

As expected, there are 45 - 15 = 30 WGBs in the trum in the 4 c 2ilR limit follows the decomposition rule 

spectrum. listed in TABLE Q] 

Af 2 (15,l,0) = -2a 2 ^, 
Af 2 (l,3,0) =4a 2 <4, 

Af 2 (6,2,+±) = 0, (B26) 
M 2 (l,l,+1) = 0, 
M 2 (1,1,0) = 4(4a + a 2 )w 2 i , 



M 2 (10,3,0) = -v 2 + 2(a + /3' 4 )uj R , 
M 2 (T0, 1, -1) = -v 2 + 2 ((a - P' A ) uj r + t)oj Ri (B27) 



It is worth noting that (1, 3, 1, 0) and (1, 1, 3, 0) remain 
degenerate which is due to the fact that for u) R — the im- 
parity is conserved by even ujbl powers. On the contrary, ^ Cj^' ^> ^) — — v + 2 (a + /34 + /? 4 ) w fl , 
in the 126# components the D-parity is broken by the r M 2 (10, 1, +1) = + 2 ((a - #4) w fl - r) oj r 
term that is linear in lobl- 



J 



M>,1,0)=( -^ + (2« + A + 2^)^ -4 72 . 2 \ 

V 2 I 2 72 ^| -z, 2 + i((4a + /3 4 + 2^)^ + 2r)^ i V 7 



The mass matrices above are defined on the bases 
{(6,l,0) s .,(6,l,0) s } and {(15, 2, -§) E . , (15, 2, 
respectively. As expected, there are in total 45 — 19 = 26 
WGBs. 



4. Few remarks on global symmetries 

It is illuminating to look at the global symmetries of 
the scalar potential when only the moduli of 45h and 
126# appear in the scalar potential. In such a case, 
(i.e., with a 2 = A 2 = A 4 = A4 = r/2 = t = /3 4 = 
/? 4 = 72 = 0), the global symmetry of Vb in Eq. ^ is 
0(45) <E> 0(252). This symmetry is spontaneously broken 
into 0(44)<g>O(251) by the 45 H and 126 ff VEVs yielding 



44 + 251 = 295 Goldstones in the scalar spectrum. Since, 
at the same time, the gauge SO(10) symmetry is broken 
to the SM gauge group, 45— 12 = 33 would-be Goldstone 
bosons are "eaten" by the gauge fields associated to the 
5O(10)/SM coset and drop from the scalar spectrum, 
295 — 33 = 262 PGB remain. Their masses are gener- 
ally expected to receive contributions from the explicitly 
breaking terms a 2 , A 2 , A 4 , A 4 , ?7 2 , t, /3 4 , (3' 4 and 72 . A 
detailed inspection of the mass matrices in Eqs. (|B5[) - 
(|B16[) . indeed, reveals the total of 262 massless degrees 
of freedom. 

In this respect, let us emphasize that the (1,3,0) and 
(8, 1, 0) components of 45# of our central interest in Sec- 
tion [TTC] belong to this category but, for various reasons, 
they receive just the a 2 -proportional mass contribution 
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at the tree level while the other couplings enter the rele- 
vant mass formulae only via loops. 



Appendix C: The tree-level gauge boson spectrum 



Let us start with the scalar kinetic tcrmtPl 



while, as expected, there are no contributions of (45#) 
to A4^(l,3,0), Mi (8, 1,0) and Mi(l,l,0). Note that, 
in the limits of the standard 5 lz (ur = ubl), flipped 
5' l z > (ur = -Ubl), 3 c 2 l 2 r 1 bl (uj r = 0) and 4 C 2 L 1 R 
(ujrl = 0) vacua, there are 25, 25, 15 and 19 massless 
gauge bosons, respectively. 



lab 



and 



1 



2-5! 



(Z? /J I])* &cde (L' AI I]) a f )c( i e , 2. Contributions to gauge bosons masses from 126h 



where the covariant derivatives are given by 



(D^ab = d^at - i-glA^ij [e^, <j)] ab 



(CI) 



(C2) 



and 



{D ^S) abcde = dpY^ a bcde i"^9i.-^-^)ij {^ij)aa'^a'bcde 

+ { e ij)bb'^ab'cde + (tij) cc'^abc' de (G3) 
~t~ {^ij^dd'^abcd'e "t - (^ij )ee' ^abcde' 



respectively, and eij 



, 10) are the 50(10) 



generators in the fundamental representation 

(<H;)ab = -i{5ai5bj - SajSbi) ■ (C4) 

One finds the following expressions for the field depen- 
dent mass matrices of the gauge bosons 

M 2 A {4>)( l3){k i) = yTr [e (y) ,<£] [t(ki),4>] (C5) 
for the contribution from the VEVs in the 45#, and 

9 2 



+ ( e (ij))bb'^ab'cde + ( e (ij) )cc' ^abc'de 
+ ( e {ij))dd'^abcd'e + ( e (ij)) ee'^abcd 

X ( e (kl))aa"^a"bcde 
+ { e (kl))bb"^ab"cde + (^(kl) )cc" ^abc"de 
+ ( e (kl))dd"^abcd"e + ( e (kl)) ee" £ abcde 



C.C. 



for that of the VEV in 126//, where (ij), (kl) stand for 
ordered pairs of indices. 



1. Contributions to gauge bosons masses from 45 h 

Evaluating Eq. (|C5[) one obtains 
^(1,1,+1)=4 S 2 ^, 

3 / — ^ ^BL ) 
(2 /<> o 5 \ „2 



M 2 A (3,l,-D^4gW B 
M 2 A (3,2,-I)=g 2 (ur-ubl) 2 
M 2 A &2,+±)=g 2 (lor + lobl? 



Besides 7^(1,3,0), X^(8,l,0) and A^(3,2,-§; 
which receive no contributions from (126#) one has 



Mi(i,i,+i) 



2 5 V| 2 

^(3,1,-|) = 2 5 2 H 2 



X 2 1 (3,2,+i) = 2 5 2 |a 




M 2 A {IXQ) 



Here the SM singlet matrix is defined on the pair of sin- 
glets from 15 and 1 of SU(4)c, respectively. One has 



DetMi(l,l,0) = 0, 
TrM A (l,l,0) = Wg 2 \a\ 



(C7) 
(C8) 



and, as expected, there is a massless state in the sin- 
glet sector. The number of vanishing entries corresponds 
to the dimension of the SU(5) algebra preserved by the 
126/f VEV a. Summing together the 45# and 126# con- 
tributions, one ends up with 12 massless vector bosons of 
the unbroken SU(3) c (^SU(2)l^U(1)y gauge symmetry. 



Appendix D: Sample spectra 

The details of the bosonic spectrum for the sample 
settings described in Sects. HTl C 1 cl and llll C 2cl are given 
in TABLES [TV] and E] respectively. The symbols RS, 
CS and VB denote, consecutively, real scalars, complex 
scalars and vector bosons. Notice that a double counting 
in the SU(2)l doublet sector is, conveniently, avoided 
by treating both (1,2, +i)'s as real scalars; for further 
details see Sect. IIII B 21 

Note also that since we are working in the Feynman 
gauge of the 5O(10) — s- SM broken phase, the Goldstones 
affect the 6-cocfficients at the mass scales of the corre- 
sponding vector bosons (recall the Feynman-gauge pole 
structure of the Goldstone boson propagators); hence, 
the Goldstone bosons (GB) are displayed along with the 
relevant gauge fields in TABLEs |IV] and [Vj 



(C6) 



19 Notice that £* 6t , df ,£ab c de = <t> a bcde<t> abcde, where £ and E* arc 
defined respectively in Eqs. (I A2I) — < T^3T > . 
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multiplet 


type 


eigenstate 


a i '-491 

A& 3 


mass [GeV] 


(6,3,+i) 


CS 


1 


(I,4,|) 


5.6 x 10 11 


(1,1,-1) 


VB 


1 


(0,0, -ii) 


1.3 x 10 14 


(1,1, +1) 


VB 


1 


(0,0, -ii) 


1.3 x 10 14 


(1,1, +1) 


GB 


1 


(0,0, i) 


1.3 x 10 14 


(1,1,0) 


VB 


1 


(0,0,0) 


2.8 x 10 14 


(1,1,0) 


GB 


1 


(0,0,0) 


2.8 x 10 14 


(8,1,0) 


RS 


1 


(1,0,0) 


7.7 x 10 14 


(3,2,+|) 


CS 


2 


(3, 2> rin) 


1.1 x 10 15 


(3,2,+!) 


CS 


1 


(3, 21 ;in) 


1.2 x 10 15 


(1,1,0) 


RS 


2 


(0,0,0) 


4.3 x 10 15 


(1,1, +2) 


CS 


1 


(o,o,|) 


4.5 x 10 15 


(3,2,-i) 


VB 


1 


( — ~y ~~, ~ 3n) 


5.2 x 10 15 


(3,2,+i) 


VB 


1 


( — T> _ T> _ 3n) 


5.2 x 10 15 


(3,2,+i) 


GB 


1 


(3' 2> 30) 


5.2 x 10 15 


(3,2, +|) 


VB 


1 


/ XI XX 55 \ 

l>— "S"' _ T> "IT J 


5.2 x 10 15 


(3,2,-|) 


VB 


1 


/ XX XX 55 \ 


5.2 x 10 15 


(3,2,-|) 


GB 


1 


/X X 5 \ 

(3' 


5.2 x 10 15 


(1,1, +1) 


CS 


2 


(0,0, i) 


5.6 x 10 15 


(1,1,0) 


RS 


3 


(0,0,0) 


5.7 x 10 15 


(1,3,0) 


RS 


1 


(0,i,0) 


6.1 x 10 15 


(3,1,+i) 


CS 


1 


G.0,4) 


6.4 x 10 15 


(8,2,+i) 


CS 


1 


(2,|,|) 


9.3 x 10 15 


(3,1,+|) 


CS 


1 


(io,i|) 


9.6 x 10 15 


(3,1,+i) 


CS 


2 


G,0,i) 


9.6 x 10 15 


(3,1,-1) 


CS 


2 


(«->o,4) 


9.6 x 10 15 


(3,1,-1) 


VB 


1 


( — — — — 1 


1.0 x 10 16 


(3,1, +2) 


VB 


1 


/ 1 1 ^ zLzL \ 

(-¥,o,-#) 


1.0 x 10 16 


(3,1,-1) 


GB 


1 


G,o,4) 


1.0 x 10 16 


(8,2,+i) 


CS 


2 


(2, it) 


1.1 x 10 16 


(6,1,+|) 


CS 


1 


(l,o,4) 


1.5 x 10 16 


(1,2,+i) 


RS 


1 


(0, T2> 2n) 


1.5 x 10 16 


(6,1,-1) 


CS 


1 


(l,o,4) 


1.5 x 10 16 


(6,1,-1) 


CS 


1 


(l,o, i) 


1.5 x 10 16 


(1.2, +3) 


RS 


2 


(0, T21 In) 


1.6 x 10 16 


(3,1,+i) 


CS 


3 


G,o,4) 


1.7 x 10 16 






1 
1 


U>2, 


l.o X 1U 


(3,2,+i) 


CS 


3 


fi i — ) 


2.1 x 10 16 


(3,2,+!) 


CS 


2 


(• 1 1 49 \ 

^ 3 > 2 > an i 


2.1 x 10 16 


(1,3,-1) 


CS 


1 


(o,l,!) 


2.6 x 10 16 


(1,1,0) 


RS 


4 


(0,0,0) 


3.0 x 10 16 



TABLE IV. A sample spectrum featuring a light (6, 3, +|) 
threshold. The relevant scalar potential parameters are given 
in the left column of Table ITTT1 A6 321 indicates the shift in the 
one-loop beta-function. The light threshold and the vector 
bosons defining the GUT scale are in boldface. 



multiplet 


type 


eigenstate 


A I '-(91 

A& 3 


mass [GeV] 


(8,2,+|) 


CS 


1 


(2,1,1) 


2.3 x 10 4 


(3,1,-|) 


VB 


1 


/ 1 1 n A A \ 

(-tt.°.-tr) 


2.8 x 10 13 


(3,1, +|) 


VB 


1 


/ 1 1 i-v AA \ 

(-¥-o,-4) 


2.8 x 10 13 


(3,1,-!) 


GB 


1 




2.8 x 10 13 


(1,1,0) 


VB 


1 


(0,0,0) 


6.1 x 10 13 


(1,1,0) 


GB 


1 


(0,0,0) 


6.1 x 10 13 


(3,2,+i) 


CS 


1 


(31 2 ' sn) 


2.6 x 10 14 


(3,2,+i) 


CS 


3 


(3, 2' sn) 


2.8 x 10 14 


(1,2,+i) 


RS 


1 


(0, 12' In) 


3.3 x 10 14 


(1,1,0) 


RS 


2 


(0,0,0) 


2.2 x 10 15 


(3,1,-1) 


CS 


2 


(*,o,4) 


2.3 x 10 15 


(6,3,+i) 


CS 


1 


(1,4,1) 


2.3 x 10 15 


(3,3,-i) 


CS 


1 


(l,2,|) 


2.3 x 10 15 


(1,3,-1) 


CS 


1 


(o,l,l) 


2.3 x 10 15 


(6,1,-|) 


CS 


1 


(l,o,4) 


3.2 x 10 15 


(1,1,0) 


RS 


3 


(0,0,0) 


3.3 x 10 15 


(8,1,0) 


RS 


1 


(1,0,0) 


4.6 x 10 15 


(1,3,0) 


RS 


1 


(0,i,0) 


6.1 x 10 15 


(3,2, +|) 


VB 


1 


/ XX XX 55 n 

\~, ~~> 


8.7 x 10 15 


(3,2,-|) 


VB 


1 


1 XX XX 55 \ 

\~, ~~, ~~> 


8.7 x 10 15 


(3,2,-|) 


GB 


1 


(a> s) 


8.7 x 10 15 


(3,2,-i) 


VB 


1 


("IT I ~3n) 


8.7 x 10 15 


(3,2,+i) 


VB 


1 


(~~, ~3n) 


8.7 x 10 15 


(3,2,+i) 


GB 


1 


(3' 2' sn) 


8.7 x 10 15 


(3,1,+i) 


CS 


1 


(l,o,4) 


1.1 x 10 16 


(3,1,+i) 


CS 


2 


(i 0,4) 


1.2 x 10 16 


(1,1, +1) 


CS 


2 


(0,0, i) 


1.6 x 10 16 


(3,1,+i) 


CS 


3 


(io,4) 


1.6 x 10 16 


(6,1,-3) 


CS 


1 


(l,o,4) 


1.6 x 10 16 


(3,2,+i) 


CS 


2 


(3' 2' 3f|) 


1.7 x 10 16 


(1,2,+i) 


RS 


2 


(0, 12' 2tl) 


1.7 x 10 16 


(8,2,+i) 


CS 


2 


(2>|,f) 


1.7 x 10 16 


(3,2,+i) 


CS 


2 


(3 ' 2 ' 30) 


1.7 x 10 16 


(1,1,-1) 


VB 


1 


(0,0, -¥) 


1.7 x 10 16 


(1,1, +1) 


VB 


1 


(0,0, -if) 


1.7 x 10 16 


[I, 1, +1) 




1 
1 


i(\ n 1 \ 
l u > u , s ) 


1.1 X 1U 


(1,1, +2) 


CS 


1 


(0,0, i) 


2.4 x 10 16 


(3,1,+|) 


CS 


1 


(- — ) 


2.4 x 10 16 


(6,1, +|) 


CS 


1 


(- — ) 


2.4 x 10 16 


(1,1,0) 


RS 


4 


(0,0,0) 


4.1 x 10 16 



TABLE V. The same as in Table HVl but for the case discussed 
in Sect. IIII C 2~cl featuring a light (8,2,+i) threshold. The 
relevant scalar potential parameters are given in the right 
column of Table fTTTl Notice that, as required by consistency, 
in both cases b S M + E Ab 321 = (-f ,-f ,-f ). 
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